Methylation-based markers for the estimation of age in African cheetah, Acinonyx jubatus.

Age is a key demographic in conservation where age classes show differences in important population metrics such as morbidity and mortality. Several traits, including reproductive potential, also show senescence with ageing. Thus, the ability to estimate age of individuals in a population is critical in understanding the current structure as well as their future fitness. Many methods exist to determine age in wildlife, with most using morphological features that show inherent variability with age. These methods require significant expertise and become less accurate in adult age classes, often the most critical groups to model. Molecular methods have been applied to measuring key population attributes, and more recently epigenetic attributes such as methylation have been explored as biomarkers for age. There are, however, several factors such as permits, sample sovereignty, and costs that may preclude the use of extant methods in a conservation context. This study explored the utility of measuring age-related changes in methylation in candidate genes using mass array technology. Novel methods are described for using gene orthologues to identify and assay regions for differential methylation. To illustrate the potential application, African cheetah was used as a case study. Correlation analyses identified six methylation sites with an age relationship, used to develop a model with sufficient predictive power for most conservation contexts. This model was more accurate than previous attempts using PCR and performed similarly to candidate gene studies in other mammal species. Mass array presents an accurate and cost-effective method for age estimation in wildlife of conservation concern.

Diversity of selected toll-like receptor genes in cheetahs (Acinonyx jubatus) and African leopards (Panthera pardus pardus).

The anthropogenic impact on wildlife is ever increasing. With shrinking habitats, wild populations are being pushed to co-exist in proximity to humans leading to an increased threat of infectious diseases. Therefore, understanding the immune system of a species is key to assess its resilience in a changing environment. The innate immune system (IIS) is the body's first line of defense against pathogens. High variability in IIS genes, like toll-like receptor (TLR) genes, appears to be associated with resistance to infectious diseases. However, few studies have investigated diversity in TLR genes in vulnerable species for conservation. Large predators are threatened globally including leopards and cheetahs, both listed as 'vulnerable' by IUCN. To examine IIS diversity in these sympatric species, we used next-generation-sequencing to compare selected TLR genes in African leopards and cheetahs. Despite differences, both species show some TLR haplotype similarity. Historic cheetahs from all subspecies exhibit greater genetic diversity than modern Southern African cheetahs. The diversity in investigated TLR genes is lower in modern Southern African cheetahs than in African leopards. Compared to historic cheetah data and other subspecies, a more recent population decline might explain the observed genetic impoverishment of TLR genes in modern Southern African cheetahs. However, this may not yet impact the health of this cheetah subspecies.

A new multiplex qPCR assay to detect and differentiate big cat species in the illegal wildlife trade.

All species of big cats, including tigers, cheetahs, leopards, lions, snow leopards, and jaguars, are protected under the Convention on the International Trade in Endangered Species (CITES). This is due in large part to population declines resulting from anthropogenic factors, especially poaching and the unregulated and illegal trade in pelts, bones, teeth and other products that are derived from these iconic species. To enhance and scale up monitoring for big cat products in this trade, we created a rapid multiplex qPCR test that can identify and differentiate DNA from tiger (Panthera tigris), cheetah (Acinonyx jubatus), leopard (Panthera pardus), lion (Panthera leo), snow leopard (Panthera uncia), and jaguar (Panthera onca) in wildlife products using melt curve analysis to identify each species by its unique melt peak temperature. Our results showed high PCR efficiency (> 90%), sensitivity (detection limit of 5 copies of DNA per PCR reaction) and specificity (no cross amplification between each of the 6 big cat species). When paired with a rapid (< 1 h) DNA extraction protocol that amplifies DNA from bone, teeth, and preserved skin, total test time is less than three hours. This test can be used as a screening method to improve our understanding of the scale and scope of the illegal trade in big cats and aid in the enforcement of international regulations that govern the trade in wildlife and wildlife products, both ultimately benefiting the conservation of these species worldwide.

A chromosome-scale high-contiguity genome assembly of the cheetah (Acinonyx jubatus).

The cheetah (Acinonyx jubatus, SCHREBER 1775) is a large felid and is considered the fastest land animal. Historically, it inhabited open grassland across Africa, the Arabian Peninsula, and southwestern Asia; however, only small and fragmented populations remain today. Here, we present a de novo genome assembly of the cheetah based on PacBio continuous long reads and Hi-C proximity ligation data. The final assembly (VMU_Ajub_asm_v1.0) has a total length of 2.38 Gb, of which 99.7% are anchored into the expected 19 chromosome-scale scaffolds. The contig and scaffold N50 values of 96.8 Mb and 144.4 Mb, respectively, a BUSCO completeness of 95.4% and a k-mer completeness of 98.4%, emphasize the high quality of the assembly. Furthermore, annotation of the assembly identified 23,622 genes and a repeat content of 40.4%. This new highly contiguous and chromosome-scale assembly will greatly benefit conservation and evolutionary genomic analyses and will be a valuable resource, e.g., to gain a detailed understanding of the function and diversity of immune response genes in felids.

Genomic analyses show extremely perilous conservation status of African and Asiatic cheetahs (Acinonyx jubatus).

We live in a world characterized by biodiversity loss and global environmental change. The extinction of large carnivores can have ramifying effects on ecosystems like an uncontrolled increase in wild herbivores, which in turn can have knock-on impacts on vegetation regeneration and communities. Cheetahs (Acinonyx jubatus) serve important ecosystem functions as apex predators; yet, they are quickly heading towards an uncertain future. Threatened by habitat loss, human-wildlife conflict and illegal trafficking, there are only approximately 7100 individuals remaining in nature. We present the most comprehensive genome-wide analysis of cheetah phylogeography and conservation genomics to date, assembling samples from nearly the entire current and past species' range. We show that their phylogeography is more complex than previously thought, and that East African cheetahs (A. j. raineyi) are genetically distinct from Southern African individuals (A. j. jubatus), warranting their recognition as a distinct subspecies. We found strong genetic differentiation between all classically recognized subspecies, thus refuting earlier findings that cheetahs show only little differentiation. The strongest differentiation was observed between the Asiatic and all the African subspecies. We detected high inbreeding in the Critically Endangered Iranian (A. j. venaticus) and North-western (A. j. hecki) subspecies, and show that overall cheetahs, along with snow leopards, have the lowest genome-wide heterozygosity of all the big cats. This further emphasizes the cheetah's perilous conservation status. Our results provide novel and important information on cheetah phylogeography that can support evidence-based conservation policy decisions to help protect this species. This is especially relevant in light of ongoing and proposed translocations across subspecies boundaries, and the increasing threats of illegal trafficking.

KIT mutations in mast cell tumours from cheetahs (Acinonyx jubatus) and domestic cats (Felis catus).

Mast cell tumours (MCT) have been documented in numerous species and mutations within the KIT proto-oncogene are implicated in the neoplastic biology of mast cells in humans, dogs and cats. This study determined high KIT gene nucleotide and Kit amino acid sequence homology between several species known to suffer mast cell neoplasia and especially high sequence conservation between the cheetah (Acinonyx jubatus) and domestic cat (Felis catus) KIT sequences. As a result, we hypothesised that KIT mutations would exist in the neoplastic DNA of four cheetahs diagnosed with MCT from a recent case series. PCR and Sanger sequencing identified conservative exon 6 KIT mutations in two of the four cheetahs. The mutations were different between the two cheetahs. Only wild-type DNA in exons 6, 8, 9 and 11 of KIT was observed in the MCTs of the remaining two cheetahs. Twenty cutaneous MCTs from domestic cats were collected for KIT mutation comparison. Twelve tumours possessed a mutation within KIT exons 6, 8 or 9 (60%, 95% CI 38.5%-81.5%). No mutations were detected in exon 11. There was no significant association between domestic feline MCT KIT mutation status and tumour histological grade (traditional schematic, P = .934; Sabattini 2-tier schematic, P = .762) or mitotic index (P = .750). KIT mRNA and Kit protein sequences are conserved across species but the role of KIT in feline MCT pathogenesis is not completely understood.

Unlocking the potential of a validated single nucleotide polymorphism array for genomic monitoring of trade in cheetahs (Acinonyx jubatus).

Cheetahs (Acinonyx jubatus) are listed as vulnerable on the International Union for Conservation of Nature Red List of Threatened Species. Threats include loss of habitat, human-wildlife conflict and illegal wildlife trade. In South Africa, the export of wild cheetah is a restricted activity under the Convention on International Trade in Endangered Species of Wild Fauna and Flora (CITES), however, limited legal trade is permitted of animals born to captive parents. To effectively monitor the legal and illegal trade in South Africa, it was thus essential to develop a validated molecular test. Here, we designed a single nucleotide polymorphism (SNP) array for cheetah from Double Digest Restriction Associated DNA sequencing data for individual identification and parentage testing. In order to validate the array, unrelated individuals and 16 family groups consisting of both parents and one to three offspring were genotyped using the Applied Biosystems™ QuantStudio™ 12K Flex Real-Time PCR System. In addition, parentage assignments were compared to microsatellite data. Cross-species amplification was tested in various felids and cheetah sub-species in order to determine the utility of the SNP array in other species. We obtained successful genotyping results for 218 SNPs in cheetah (A. j. jubatus) with an optimal DNA input concentration ranging from 10 to 30 ng/µl. The combination of SNPs had a higher resolving power for individual identification compared to microsatellites and provided high assignment accuracy in known pedigrees. Cross-species amplification in other felids was determined to be limited. However, the SNP array demonstrated a clear genetic discrimination of two cheetah subspecies tested here. We conclude that the described SNP array is suitable for accurate parentage assignment and provides an important traceability tool for forensic investigations of cheetah trade.

Ancient mtDNA from the extinct Indian cheetah supports unexpectedly deep divergence from African cheetahs.

The Indian cheetah was hunted to extinction by the mid-20th century. While analysis of 139 bp of mitochondrial DNA (mtDNA) has confirmed that the Indian cheetah was part of the Asiatic subspecies (Acinonyx jubatus venaticus), the detailed relationships between cheetah populations remains unclear due to limited genetic data. We clarify these relationships by studying larger fragments of cheetah mtDNA, both from an Indian cheetah museum specimen and two African cheetah, one modern and one historic, imported into India at different times. Our results suggest that the most recent common ancestor of cheetah mtDNA is approximately twice as ancient as currently recognised. The Indian and Southeast African (Acinonyx jubatus jubatus) cheetah mtDNA diverged approximately 72 kya, while the Southeast and Northeast African (Acinonyx jubatus soemmeringii) cheetah mtDNA diverged around 139 kya. Additionally, the historic African cheetah sampled from India proved to have an A. j. jubatus haplotype, suggesting a hitherto unrecognised South African route of cheetah importation into India in the 19th century. Together, our results provide a deeper understanding of the relationships between cheetah subspecies, and have important implications for the conservation of A. j. venaticus and potential reintroduction of cheetahs into India.

Conservation Genetics of the Cheetah: Lessons Learned and New Opportunities.

The dwindling wildlife species of our planet have become a cause célèbre for conservation groups, governments, and concerned citizens throughout the world. The application of powerful new genetic technologies to surviving populations of threatened mammals has revolutionized our ability to recognize hidden perils that afflict them. We have learned new lessons of survival, adaptation, and evolution from viewing the natural history of genomes in hundreds of detailed studies. A single case history of one species, the African cheetah, Acinonyx jubatus, is here reviewed to reveal a long-term story of conservation challenges and action informed by genetic discoveries and insights. A synthesis of 3 decades of data, interpretation, and controversy, capped by whole genome sequence analysis of cheetahs, provides a compelling tale of conservation relevance and action to protect this species and other threatened wildlife.

The difference between trivial and scientific names: There were never any true cheetahs in North America.

Dobrynin et al. (Genome Biol 16:277, 2015) recently published the complete genome of the cheetah (Acinonyx jubatus) and provided an exhaustive set of analyses supporting the famously low genetic variation in the species, known for several decades. Their genetic analyses represent state-of-the-art and we do not criticize them. However, their interpretation of the results is inconsistent with current knowledge of cheetah evolution. Dobrynin et al. suggest that the causes of the two inferred bottlenecks at ∼ 100,000 and 10,000 years ago were immigration by cheetahs from North America and end-Pleistocene megafauna extinction, respectively, but the first explanation is impossible and the second implausible.

Serum Amyloid A Protein Concentration in Blood is Influenced by Genetic Differences in the Cheetah (Acinonyx jubatus).

Systemic amyloid A (AA) amyloidosis is a major cause of morbidity and mortality among captive cheetahs. The self-aggregating AA protein responsible for this disease is a byproduct of serum amyloid A (SAA) protein degradation. Transcriptional induction of the SAA1 gene is dependent on both C/EBPß and NF-κB cis-acting elements within the promoter region. In cheetahs, 2 alleles exist for a single guanine nucleotide deletion in the putative NF-κB binding site. In this study, a novel genotyping assay was developed to screen for the alleles. The results show that the SAA1A (-97delG) allele is associated with decreased SAA protein concentrations in the serum of captive cheetahs (n = 58), suggesting genetic differences at this locus may be affecting AA amyloidosis prevalence. However, there was no significant difference in the frequency of the SAA1A (-97delG) allele between individuals confirmed AA amyloidosis positive versus AA amyloidosis negative at the time of necropsy (n = 48). Thus, even though there is evidence that having more copies of the SAA1A (-97delG) allele results in a potentially protective decrease in serum concentrations of SAA protein in captive cheetahs, genotype is not associated with this disease within the North American population. These results suggest that other factors are playing a more significant role in the pathogenesis of AA amyloidosis among captive cheetahs.

Genomic legacy of the African cheetah, Acinonyx jubatus.

BACKGROUND: Patterns of genetic and genomic variance are informative in inferring population history for human, model species and endangered populations. RESULTS: Here the genome sequence of wild-born African cheetahs reveals extreme genomic depletion in SNV incidence, SNV density, SNVs of coding genes, MHC class I and II genes, and mitochondrial DNA SNVs. Cheetah genomes are on average 95 % homozygous compared to the genomes of the outbred domestic cat (24.08 % homozygous), Virunga Mountain Gorilla (78.12 %), inbred Abyssinian cat (62.63 %), Tasmanian devil, domestic dog and other mammalian species. Demographic estimators impute two ancestral population bottlenecks: one >100,000 years ago coincident with cheetah migrations out of the Americas and into Eurasia and Africa, and a second 11,084-12,589 years ago in Africa coincident with late Pleistocene large mammal extinctions. MHC class I gene loss and dramatic reduction in functional diversity of MHC genes would explain why cheetahs ablate skin graft rejection among unrelated individuals. Significant excess of non-synonymous mutations in AKAP4 (p<0.02), a gene mediating spermatozoon development, indicates cheetah fixation of five function-damaging amino acid variants distinct from AKAP4 homologues of other Felidae or mammals; AKAP4 dysfunction may cause the cheetah's extremely high (>80 %) pleiomorphic sperm. CONCLUSIONS: The study provides an unprecedented genomic perspective for the rare cheetah, with potential relevance to the species' natural history, physiological adaptations and unique reproductive disposition.

Development of MHC-Linked Microsatellite Markers in the Domestic Cat and Their Use to Evaluate MHC Diversity in Domestic Cats, Cheetahs, and Gir Lions.

Diversity within the major histocompatibility complex (MHC) reflects the immunological fitness of a population. MHC-linked microsatellite markers provide a simple and an inexpensive method for studying MHC diversity in large-scale studies. We have developed 6 MHC-linked microsatellite markers in the domestic cat and used these, in conjunction with 5 neutral microsatellites, to assess MHC diversity in domestic mixed breed (n = 129) and purebred Burmese (n = 61) cat populations in Australia. The MHC of outbred Australian cats is polymorphic (average allelic richness = 8.52), whereas the Burmese population has significantly lower MHC diversity (average allelic richness = 6.81; P < 0.01). The MHC-linked microsatellites along with MHC cloning and sequencing demonstrated moderate MHC diversity in cheetahs (n = 13) and extremely low diversity in Gir lions (n = 13). Our MHC-linked microsatellite markers have potential future use in diversity and disease studies in other populations and breeds of cats as well as in wild felid species.

Immunogenetic variation and differential pathogen exposure in free-ranging cheetahs across Namibian farmlands.

BACKGROUND: Genes under selection provide ecologically important information useful for conservation issues. Major histocompatibility complex (MHC) class I and II genes are essential for the immune defence against pathogens from intracellular (e.g. viruses) and extracellular (e.g. helminths) origins, respectively. Serosurvey studies in Namibian cheetahs (Acinonyx juabuts) revealed higher exposure to viral pathogens in individuals from north-central than east-central regions. Here we examined whether the observed differences in exposure to viruses influence the patterns of genetic variation and differentiation at MHC loci in 88 free-ranging Namibian cheetahs. METHODOLOGY/PRINCIPAL FINDINGS: Genetic variation at MHC I and II loci was assessed through single-stranded conformation polymorphism (SSCP) analysis and sequencing. While the overall allelic diversity did not differ, we observed a high genetic differentiation at MHC class I loci between cheetahs from north-central and east-central Namibia. No such differentiation in MHC class II and neutral markers were found. CONCLUSIONS/SIGNIFICANCE: Our results suggest that MHC class I variation mirrors the variation in selection pressure imposed by viruses in free-ranging cheetahs across Namibian farmland. This is of high significance for future management and conservation programs of this species.

Specifying and sustaining pigmentation patterns in domestic and wild cats.

Color markings among felid species display both a remarkable diversity and a common underlying periodicity. A similar range of patterns in domestic cats suggests a conserved mechanism whose appearance can be altered by selection. We identified the gene responsible for tabby pattern variation in domestic cats as Transmembrane aminopeptidase Q (Taqpep), which encodes a membrane-bound metalloprotease. Analyzing 31 other felid species, we identified Taqpep as the cause of the rare king cheetah phenotype, in which spots coalesce into blotches and stripes. Histologic, genomic expression, and transgenic mouse studies indicate that paracrine expression of Endothelin3 (Edn3) coordinates localized color differences. We propose a two-stage model in which Taqpep helps to establish a periodic pre-pattern during skin development that is later implemented by differential expression of Edn3.

The clavicles of Smilodon fatalis and Panthera atrox (Mammalia: Felidae) from Rancho La Brea, Los Angeles, California.

The Rancho La Brea collections at the George C. Page Museum in Los Angeles, California, contain the largest single inventory of Smilodon fatalis remains representing virtually every bone in the skeleton. Eighteen clavicles of two distinctive shapes have been recovered from historical and recent excavations at Rancho La Brea. In this study, we identify these specimens to species through comparison of their morphology and morphological variability with clavicles found in modern felids. This study includes a reevaluation of clavicles that have been previously assigned to S. fatalis, which are more likely to be those of Panthera atrox, and the description of pantherine cat clavicles. A previously undescribed sample of clavicles not only includes some of the same pantherine morph but also 10 specimens, herein assigned to S. fatalis, which are morphologically distinctive and significantly smaller than the previously described specimens. In addition, we report unexpected variations between clavicles of Panthera leo and P. tigris: the clavicles of P. leo closely resemble those of the large Rancho La Brea clavicle morph-which presumably belongs to P. atrox-thus supporting a P. leo/P. atrox clade. We report distinctive morphology of the clavicles of Acinonyx jubatus. Possible functional and phylogenic significance of felid clavicles is suggested.

Cheetahs have 4 serum amyloid a genes evolved through repeated duplication events.

Amyloid A (AA) amyloidosis is a leading cause of mortality in captive cheetahs (Acinonyx jubatus). We performed genome walking and PCR cloning and revealed that cheetahs have 4 SAA genes (provisionally named SAA1A, SAA1B, SAA3A, and SAA3B). In addition, we identified multiple nucleotide polymorphisms in the 4 SAA genes by screening 51 cheetahs. The polymorphisms defined 4, 7, 6, and 4 alleles for SAA1A, SAA3A, SAA1B, and SAA3B, respectively. Pedigree analysis of the inheritance of genotypes for the SAA genes revealed that specific combinations of alleles for the 4 SAA genes cosegregated as a unit (haplotype) in pedigrees, indicating that the 4 genes were linked on the same chromosome. Notably, cheetah SAA1A and SAA1B were highly homologous in their nucleotide sequences. Likewise, SAA3A and SAA3B genes were homologous. These observations suggested a model for the evolution of the 4 SAA genes in cheetahs in which duplication of an ancestral SAA gene first gave rise to SAA1 and SAA3. Subsequently, each gene duplicated one more time, uniquely making 4 genes in the cheetah genome. The monomorphism of the cheetah SAA1A protein might be one of the factors responsible for the high incidence of AA amyloidosis in this species.

Digital gene expression for non-model organisms.

Next-generation sequencing technologies offer new approaches for global measurements of gene expression but are mostly limited to organisms for which a high-quality assembled reference genome sequence is available. We present a method for gene expression profiling called EDGE, or EcoP15I-tagged Digital Gene Expression, based on ultra-high-throughput sequencing of 27-bp cDNA fragments that uniquely tag the corresponding gene, thereby allowing direct quantification of transcript abundance. We show that EDGE is capable of assaying for expression in >99% of genes in the genome and achieves saturation after 6-8 million reads. EDGE exhibits very little technical noise, reveals a large (10(6)) dynamic range of gene expression, and is particularly suited for quantification of transcript abundance in non-model organisms where a high-quality annotated genome is not available. In a direct comparison with RNA-seq, both methods provide similar assessments of relative transcript abundance, but EDGE does better at detecting gene expression differences for poorly expressed genes and does not exhibit transcript length bias. Applying EDGE to laboratory mice, we show that a loss-of-function mutation in the melanocortin 1 receptor (Mc1r), recognized as a Mendelian determinant of yellow hair color in many different mammals, also causes reduced expression of genes involved in the interferon response. To illustrate the application of EDGE to a non-model organism, we examine skin biopsy samples from a cheetah (Acinonyx jubatus) and identify genes likely to control differences in the color of spotted versus non-spotted regions.

Phylogeography, genetic structure and population divergence time of cheetahs in Africa and Asia: evidence for long-term geographic isolates.

The cheetah (Acinonyx jubatus) has been described as a species with low levels of genetic variation. This has been suggested to be the consequence of a demographic bottleneck 10 000-12 000 years ago (ya) and also led to the assumption that only small genetic differences exist between the described subspecies. However, analysing mitochondrial DNA and microsatellites in cheetah samples from most of the historic range of the species we found relatively deep phylogeographic breaks between some of the investigated populations, and most of the methods assessed divergence time estimates predating the postulated bottleneck. Mitochondrial DNA monophyly and overall levels of genetic differentiation support the distinctiveness of Northern-East African cheetahs (Acinonyx jubatus soemmeringii). Moreover, combining archaeozoological and contemporary samples, we show that Asiatic cheetahs (Acinonyx jubatus venaticus) are unambiguously separated from African subspecies. Divergence time estimates from mitochondrial and nuclear data place the split between Asiatic and Southern African cheetahs (Acinonyx jubatus jubatus) at 32 000-67 000 ya using an average mammalian microsatellite mutation rate and at 4700-44 000 ya employing human microsatellite mutation rates. Cheetahs are vulnerable to extinction globally and critically endangered in their Asiatic range, where the last 70-110 individuals survive only in Iran. We demonstrate that these extant Iranian cheetahs are an autochthonous monophyletic population and the last representatives of the Asiatic subspecies A. j. venaticus. We advocate that conservation strategies should consider the uncovered independent evolutionary histories of Asiatic and African cheetahs, as well as among some African subspecies. This would facilitate the dual conservation priorities of maintaining locally adapted ecotypes and genetic diversity.